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Spin-based electronics or spintronics relies on the ability to store, transport and manipulate electron 
spin polarization with great precision [D 121 [3], [4j. In its ultimate limit, information is stored in the 
spin state of a single electron, at which point also quantum information processing becomes a possibil¬ 
ity HHU- Here we demonstrate the manipulation, transport and read-out of individual electron spins 
in a linear array of three semiconductor quantum dots. First, we demonstrate single-shot read-out 
of three spins with fidelities of 97% on average, using an approach analogous to the operation of a 
charge-coupled-device (CCD) [3- Next, we perform site-selective control of the three spins thereby 
writing the content of each pixel of this “Single-Spin CCD”. Finally, we show that shuttling an electron 
back and forth in the array hundreds of times, covering a cumulative distance of 80 /im, has negligible 
influence on its spin projection. Extrapolating these results to the case of much larger arrays, points 
at a diverse range of potential applications, from quantum information to imaging and sensing. 

Past experiments have shown macroscopic spin transport over distances exceeding 100 pm in clean bulk 2D semi¬ 
conductors HUS]- Furthermore, controlled single electron transport through semiconductor nanostructures m is now 
routine, with applications ranging from current standards to sensors and digital electronics DU- Also the controlled 
transfer of individual electrons between nanostructures separated by several micron has been realized [1211131 . However, 
the key combination of single-electron transport and spin preservation over large distances remains to be demonstrated. 

A promising platform for moving around individual spins in a controlled manner is provided by analogy to a 
CCD |7J. In a CCD, pockets of electrical charge are passed on along a capacitor array that acts as a shift-register, 
similar in spirit to a bucket-brigade. The pockets of charge arrive sequentially at the end of the array, where they 
are detected via a charge amplifier. This simple concept has enabled CCD cameras containing millions of pixels with 
applications from consumer electronics to astronomy jl4j . The creation of an analogous device that can operate and 
detect single spins would have powerful and diverse applications as well. For instance, it would not only enable reading 
out the outcome of a large quantum simulation or computation performed in a 2D array of spins mm, but also could 
be used for coherent imaging at the single photon level HE] or magnetic field sensing with 200 nm spatial resolution, 
using each single spin as a local probe. We term such a device a Single-Spin CCD. Its operation requires the ability to 
shuttle spins one by one along a chain of sites without disturbing the spin states, and to record the state of the spin 
at the end of the chain. 

We have created a prototype Single-Spin CCD using a linear triple quantum dot array. The array is formed electro¬ 
statically in a 2D electron gas (2DEG) 85 nm below the surface, see Fig. la. Gate electrodes fabricated on the surface 
of a GaAs/AlGaAs heterostructure (see Methods) are biased with appropriate voltages to selectively deplete regions 
of the 2DEG and define the linear array of three quantum dots. Each dot is initially occupied with one electron, and 
each of the three electrons constitutes a single spin-1 particle. The main function of gates LS and RS is to set the 
tunnel coupling with the left and right reservoir, respectively. D\ and D 2 control the interdot tunnel coupling (tuned 
to 0.8 and 0.5 GHz respectively) and Pi, P 2 and P 3 are used to set the electron number in each dot. The sensing dot 
(SD) next to the quantum dot array is used for non-invasive charge sensing using radiofrequency (RF) reflectometry 
to monitor the number of electrons in each dot |18j . An in-plane magnetic field B ext = 3.51 T is applied to split the 
spin-up (t) and spin-down (j.) states of each electron by the Zeeman energy (E z ~ 87 ^eV) defining a qubit in each 
of the dots. The electron temperature is 75 mK. 
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The CCD is initialized by loading f-spins from the right reservoir to every dot as described by the pulse sequence 
depicted by the blue and red arrows in Figs. lbc. Once the desired spin manipulation has been performed (see below) in 
the (l,l,l)-regime, the read-out sequence is started. This sequence follows the reverse path from the loading sequence, 
with the addition of three read-out positions that are denoted by green circles. As in a CCD, the three electrons are 
pushed sequentially to the read-out site at the end of the chain. First, the right dot is tuned to the position of green 
circle nr. 3 in Fig. lc. At this position, an excited spin-), electron is allowed to tunnel to the reservoir, whilst a ground 
state spin-f electron will remain in the dot. The nearby sensing dot is used to record whether or not the electron 
tunneled out, revealing its spin state m■ Next, we adjust the gate voltages to shuttle the center electron to the right 
dot for read out, at the position of green circle nr. 2. Afterwards we shuttle the left spin through the center dot to 
the right, and complete the three-spin read-out at the position of green circle nr. 1 in Fig. lb. Details on the pulse 
scheme can be found in the Supplementary Information. 

We test the operation of the Single-Spin CCD by preparing various combinations of the eight three-spin popula¬ 
tions ftt through HI. To do so, we implement site-selective manipulation of each of the three spins exploiting the 
small difference in ^-factors between the dots (Fig. 2i). We chirp a microwave electric held from well below to well 
above one of the three spin resonance frequencies, which results in adiabatic inversion of the corresponding spin via 
electric dipole spin resonance (EDSR) [5], here with 76% spin-), initialization efficiency on average (see Supplementary 
Information). This amounts to ‘writing’ the qubits of the Single-Spin CCD. Starting from ftti we create in this way 
the spin states j,ft> tit or ttl and subsequently vary the waiting time in the (1,1,^-configuration. During that time, 
the populations evolve as the spins relax back to the ground state ttt- 

Figs. 2a-h show the eight three-spin-probabilities as a function of the waiting time for the three initial state prepa¬ 
rations. In addition, data are shown where random spins are injected in each of the three dots. We see that the data 
follow closely the expected behavior (see caption), shown as solid lines, indicating proper operation of the Single-Spin 
CCD. More specifically, within the bounds of our measurement fidelities and statistical errors, there is no evidence of 
one measurement influencing the next, nor of mixing between adjacent spins driven by the exchange interaction [3j. 
Such mixing would be visible as an initial rise in the .(.-fraction for a spin when its neighbor was prepared as I (for 
instance, the |tt and ftl populations would initially rise for the case where fit is prepared). 

This is the first demonstration of reading out multiple spins through the same reservoir. Read-out fidelities are on 
average 98.2 (±0.5)% for spin-up and 95.8 (±0.3)% for spin-down (see Supplementary Information for details). The 
spin-down fidelities are mainly limited by the Tj-decay (~10 ms) of the spins waiting to be read out and the finite 
measurement bandwidth through which spin-), tunnel events are sometimes missed. The spin-up fidelities are likely 
limited by thermal broadening and/or electrical noise. 

Next we examine the effect of shuttling electrons between dots on their spin projection. We anticipate three 
mechanisms that could in principle cause the spin-projection to change: (i) charge exchange with the reservoirs, (ii) 
spin-orbit (SO) interaction and (iii) hyperfine interaction with the nuclear spins of the quantum dot host material. 
Exchange with the reservoirs (i) is suppressed by applying precisely tuned pulse sequences and keeping the tunnel 
barriers of the reservoirs sufficiently closed. For the present gate voltage settings, we estimate this error to be < 10~ 5 
per hop along the array (see Supplementary Information). The SO-interaction (ii) could affect the spin state, in a 
deterministic way, as it propagates through the array [22l [23] . The direction of movement with respect to the crystal 
axis determines the magnitude and direction of the SO-field. We expect the SO-interaction to be largest for motion 
along the [110] axis and minimal along the interdot axis [110] |24lI25| . see (Fig. la). Furthermore, the small SO-field 
still originating from movement along this interdot axis will be parallel to the external field and therefore will not 
influence the spin projection. We thus expect the effect of SO to be very small. The hyperfine interaction (iii) can 
cause random flips arising from the instantaneous unknown difference in perpendicular (relative to B ext ) hyperfine 
field, 5B±, between neighbouring dots. This effect is suppressed by the large B ext and estimated to be < 1 • 10 -6 per 
shuttle event assuming SB± < 7 mT (see Supplementary Information), and could even be corrected for using real-time 
Hamiltonian estimation [26] . 

To verify experimentally whether it is possible to shuttle electrons while preserving their spin projection, we simu¬ 
late a very large array using the triple dot device. Using the charge states from Fig. lb, we load one random electron 
in the (0,0,l)-state and shuttle it back and forth many times to (1,0,0) by passing through (0,1,0) as depicted schemat¬ 
ically in Fig. 3b. Each run of going back and forth constitutes a total of four jumps from one dot to the other. The 
time between hops is kept much longer than the single-spin Tff, which prevents coherent error accumulation in phase 
space. The effect of moving back and forth is therefore the same as traversing four dots in the same direction. In 
Fig. 3a, we vary both the total time spent in the CCD, tccD , and the number of interdot hops, rihops-, and read out 
the spin at the end. Fitting the data for each tccD to a linear curve gives an average change in the spin-down fraction 
ranging between —1.7 • 10 -6 and ±8.3 • 10~ 6 per hop (see Supplementary Information). We compare the measured 
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spin state after shuttling many times, with the spin state measured after the electron has shuttled back and forth only 
once in exactly the same total amount of time. The latter effectively constitutes a weighted T\ measurement over the 
three dots. The fact that the triangular symbols in Fig. 3a fall exactly on the weighted Xj decay, indicates that there 
is no sign of spin flips other than through spontaneous relaxation, even after more than 500 hops. Taking an interdot 
distance of 160 nm, this corresponds to a total distance of about 80 /im. 

An important question to address is the scalability of this approach. The data of Fig. 3a shows that in the current 
experiment, the bottleneck is not the shuttling of the spins itself, but rather relaxation in the time the spins have to 
wait before they can be read out. More specifically, the limiting factors for this experiment are the times it takes to 
read out (~130/zs), empty a dot (^75p,s) and the time allowed for shuttling to the next dot dot (~10/is). The black 
curves in Fig. 3c show the predicted spin-down fidelity as a function of the CCD-length, extrapolated based on the 
present numbers. 

With a few technical improvements involving additional pulse lines, larger interdot tunnel couplings, and a some¬ 
what lower magnetic field, tunneling and emptying can occur on a ns timescale, the read-out time can be halved and 
the Ti can be doubled (see Supplementary Information). This would give fidelities as shown in red in Fig. 3c. They 
allow one to read out 50 spins with > 83% fidelity. Adding another SD on the other side of the array will double the 
possible CCD length. 

The read-out time can be shortened more dramatically via the inclusion of Pauli spin blockade (PSB) in the read¬ 
out scheme. PSB allows one to distinguish whether two neighboring spins are parallel or anti-parallel E3 HB]. psb 
read-out within 1 /is for a fidelity of 97% has been demonstrated in GaAs dots [25] . To implement this method we use 
the right spin as an ancillary spin and one additional dot to quickly initialize this spin in the f-state using a so-called 
hot-spot where spins relax on a sub-fxs timescale ISSUE]. Then we compare spin 1 and 2 using PSB and from that 
determine the state of spin 2. After discarding spin 1, we initialize spin 2 to spin-up, again using a hot-spot, and then 
bring it close to the charge detector by shuttling it to the rightmost dot. Spin 3 can now occupy dot 2 and its spin-state 
can be determined using PSB etc. This leads to the blue curve in Fig. 3c predicting fidelities above 88 % for arrays 
larger than 1000 spins. Moving to a different host material such as Si or SiGe, Tf-times can reach seconds mm, 
boosting fidelities further (green curve). For such large Tj-times, we estimate that the read-out fidelity of the last spin 
that is read out in an array of 1000 dots, will decrease by only 0.07% compared to the read-out fidelity of the first 
spin. 

The high fidelity with which the spin projection can be preserved upon shuttling between dots thus allows scal¬ 
ing the Single-Spin CCD concept to linear arrays of hundreds of dots. Two-dimensional arrays would be a natural 
next step. Besides innovative gate designs, this would either require the use of weak spin-orbit materials like silicon, 
or (predictable) spin-orbit induced rotations need to be taken into account along at least one direction. Finally, in 
materials with negligible hyperfine coupling, such as 28 Si enriched substrates [33], not only the spin projection but 
also the spin phase is expected to be preserved during shuttling. Such coherent single-spin shuttles allow qubits to be 
moved in the course of a quantum computation, an essential tool for powerful quantum computing architectures |3SE]- 
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Figure 1 Linear array of three quantum dots and Single-Spin CCD operation, a SEM image of a sample 
nominally identical to the one used for the measurements. Dotted circles indicate quantum dots, squares indicate Fermi 
reservoirs in the 2DEG, which are connected to ohmic contacts. The RF reflectance of the SD is monitored in order to 
determine the occupancies of the triple dot. b,c Charge stability diagram of the triple dot for two different values of 
M (b M = -42 mV, c M = -56 mV). L, M and R are linear combinations of the voltages applied to gates Pi, P 2 and 
P3, allowing us to partially compensate for cross-capacitances (see Supplementary Information). The occupancy of 
each dot is denoted by (n, m,p) corresponding to the number of electrons in the left, middle and right dot respectively. 
The fading of the middle dot charge transition lines can be explained in a similar way as in |85j (black dashed lines 
indicate their positions). The CCD is initialized by loading three electrons from the right reservoir following first the 
blue arrows in b, and then pulsing M to continue loading along the red arrows in c (analogous to a shift-register). 
We load t-spins by ramping slowly through the charge transition lines with the right reservoir (slow compared to the 
tunnel rate with the reservoir). Read-out occurs at the position of the green circles using spin-selective tunneling (see 
inset of c). 
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Figure 2 Writing and reading out the Single-Spin CCD. a-h, Measured three-spin populations as a function 
of waiting time between state preparation and read-out for four different state preparations. Starting from ftt (> 95% 
initialization efficiency), states 4-tt and ttl are prepared by EDSR in the (l,l,l)-regime. State f4-t is created by 
applying EDSR in the (l,l,0)-regime, and then loading the third f-spin. For the fourth state preparation, the gate 
voltages are not ramped but pulsed across the relevant charge transitions, so that electrons with random spin will 
occupy the three dots. The contrast is limited to 0.8 mostly by imperfect adiabatic inversion and not by the read-out 
fidelities. Each datapoint is an average of 999 measurements (error bars two s.d.). Solid lines are products of the 
calculated single-spin probabilities based on the individual Tf and initial spin-down probability of each dot. The fact 
that they overlap with the corresponding three-spin probabilities demonstrates that no (unintended) correlations were 
introduced between the spins, i Resonance frequency for each dot * as a function of the magnetic field. A fit of 
the form f res = 9l ^ B -f [7], with the Bohr magneton, h Plancks constant, gives g\ = —0.430 ± 0.001, 

g\ = —0.434 ± 0.003 and g\ = —0.434 ± 0.002. Despite similar (q-factors in dot 2 and 3, time-variations of the local 
nuclear field still give rise to stable configurations where we can selectively address the two dots (see Supplementary 
Information, also for g\). 
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Figure 3 Preservation of the spin-projection during shuttling, a Circles and squares: measured spin- 
down probability after nhops interdot hopping events for three values of the total time tcco ■ Black dashed lines 
are linear fits to the data. Diamonds: measured spin-down probability after shuttling back and forth just once, as 
a function of the total shuttling time, tcGD- The decay time constant represents a weighted Tj over all three dots 
weighted = 10-6 ± 0.5 ms). Triangles depict the average value of the shuttling data (circles and squares) for each 
value of tccD (error bars two s.d.). b Schematic representation of the tunneling of a spin back and forth inside the 
CCD array. Arrows depict the backwards trajectory from the right to the left dot. Reversing this pathway returns the 
spin to the right dot. There are two separate stages with the electron in the middle dot, to prevent charge exchange 
with the reservoirs (see Supplementary Information), c Estimated spin-down measurement fidelity as a function of 
the length of the Single-Spin CCD for the current settings, improved conditions as described in the main text and for 
the proposed PSB-scheme in GaAs and SiGe. Solid lines indicate the fidelity for the spin that is read out last, dashed 
lines the average fidelity for the whole array. The spin-up fidelity is independent of the CCD size. 
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1 Methods and materials 


The experiment was performed on a GaAsZAlo. 25 Gao. 75 As heterostructure grown by molecular-beam epitaxy, with 
a 85-nm-deep 2DEG with an electron density of 2.0 • 10 11 cm -2 and mobility of 5.6 • 10 6 cm 2 V” 1 s _1 at 4 K. The 
metallic (TiAu) surface gates were fabricated using electron-beam lithography. The device was cooled inside an 
Oxford Kelvinox 400HA dilution refrigerator to a base temperature of 45 mK. To reduce charge noise the sample was 
cooled while applying a positive voltage on all gates (ranging between 0 and 350 mV) [Tj. The device was tuned in 
the single-electron regime. The tunnel coupling at zero detuning between dot 1 and 2 was measured to be 0.8 GHz 
and 0.5 GHz between dot 2 and 3 using photon-assisted tunneling [2|. Gates Pl, P 2 , P 3 and D 2 were connected to 
homebuilt bias-tees (i?C=470 ms), enabling application of d.c. voltage bias as well as high-frequency voltage excitation 
to these gates. RF reflectometry of the SD was performed using an LC circuit matching a carrier wave of frequency 
110.35 MHz. The inductor is formed by a microfabricated NbTiN superconducting spiral inductor with an inductance 
of 3.0/iH. The power of the carrier wave arriving at the sample was estimated to be -99 dBm. The carrier signal 
is only unblanked during read-out. The reflected signal was amplified using a cryogenic Weinreb CITLF2 amplifier 
and subsequently demodulated using homebuilt electronics. Real time data acquisition was performed using a FPGA 
DEO-Nano Terasic programmed to detect tunnel events using a Schmidt trigger. The microwaves were generated using 
a HP83650A connected to P 2 via a homemade bias-tee at room temperature. Voltage pulses to the gates were applied 
using a Tektronix AWG5014 (0-100% risetime = 5 ns). 

2 Calculation of the fidelities 

In this part we discuss how we obtain the read-out fidelities from the data. We follow a similar approach as discussed 
in [3] but extend the method to more dots, and incorporate some minor differences to comply with the CCD-scheme. 
We denote single spin read-out fidelity as F- where j denotes the dot in order of the read-out (1 is the rightmost dot 
for the current sample, 2 is the center dot, and 3 the leftmost dot), and i £ {{,{}■ 

2.1 Analytic expressions for the fidelity 

We start with the fidelity for the spin-up state, Ff. = 1 — a?, where a? is the probability that a step is detected in 
the SD signal although there was a spin-up electron in the quantum dot. This can occur if the SD signal exceeds 
the threshold even though the electron stayed in the dot (referred to as a signal processing error), or if the electron 
tunnels out of the dot due to thermal or electric field fluctuations. We can determine a 3 directly for each quantum 
dot by initializing in the spin-up state and successively reading it out. In this case, the fraction of the cases where 
a spin-up electron is declared spin-down directly gives a 3 . This method assumes perfect initialization and thus gives 
an upper bound on a?. op might be slightly different for the different dots due to differences in the tunnel barrier 
with the right reservoir for each charge state, the length of the read-out stages, the thresholds and the signal heights. 
Except for these slight differences, Fi is independent of the size of the array. There are two effects that may influence 
a spin-up state during the read-out sequence: (1) back-action from the SD or (2) thermal excitation. In the CCD 
scheme, each electron will experience a similar amount of back-action from the SD (an electron far away in an array 
will not yet experience back-action). Thermal excitation is negligible for the current energy scales: E z ss 90 /reV and 
ksT ss 6.5 [ieV. In thermal equilibrium, an initially perfect spin-up electron will have a probability of ~ 1CV 6 to be 
in the spin-down state (Boltzmann distribution). This is a negligible correction for our current spin-up fidelities. 

Next, we estimate the fidelity for the spin-down state based on the relaxation time T(, the duration of the read-out 
stage, T^, the electron tunneling rates (denoted by T 1 ^, where a £ {in, out} represents tunneling in or out of the 
quantum dot) and an analysis of the measurement bandwidth. 

The spin-down fidelity for the first dot, j = 1, is if 1 = S 1 + A 1 . 'EP describes the probability that a spin-down 
electron tunnels out during the read-out stage before it relaxes. hP describes the probability that a spin-down electron 
relaxes during the read-out stage, but still tunnels out because of the nonzero T. 

First we assume infinite measurement bandwidth. In that case E J can be defined as follows: 

t h 

V = J r iMt)dt (1) 

0 

where T^P^t) is the probability density function for a spin-down electron to tunnel out at time t. The proba- 

dP 

bility for the electron to be spin-down, T}(t), or spin-up, P-\{t), at time t follows from rate equations: (1) = 

—T^Pf + (l/Ti) P^ and (2) ^ = — (T^ t + l/Ti)P±. These equations describe the probabilities before the read-out 
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event occurred (after which a spin-up electron will fill the dot with rate rfrf ). The initial conditions are P^(t = 0) = 0 
and Prft = 0) = 1. From this we find P-j-(t) = — d —(e _r °«A — e _ ^ T i +r °"‘^) and P^(f) = e _ ^ T i +r °“ t ^. 

Next we add the measurement bandwidth to the calculation. A tunnel event leads to a pulse in the SD signal. 
The duration of this pulse is determined by T rf . A pulse with duration r will be detected with probability Brfr), 
which depends on the measurement bandwidth (see below). Therefore the probability that a pulse which is caused by 
an electron that tunnels out is detected is: 


T - 7 

1 R OO 

= / I rfrfPrft)dt j Trfe^^Brfrfdr (2) 

o o 

We can extend the limit of the second integral to oo because we also count events where an electron did not yet 
tunnel back during the read-out time as spin-down. Eq. Q assumes that we detect every electron leaving the dot with 
a 100% probability. This assumption is not correct for tunnel events occurring very close to the end of the read-out 
stage. Based on the measured bandwidth (see below) we can assume that every electron leaving the dot at least 0.7 
/is before the end of the read-out stage will be detected. This gives the following lower bound for S- 7 : 

T^ —0.7 fis oo 

= J I T^e- r ^ T B\r)dr (3) 

o o 

At the end of each read-out we always empty the right dot before going to the next read-out position so subsequent 
read-out stages always have the same charge configuration. 

The probability that a spin-down electron relaxes during the read-out stage before it could tunnel out but that 
nevertheless a step is detected due to the spin-up electron tunneling out is given by 

Tr— 0-7/is oo 

A 7 = J r rf t Prft)dt J T^e- r ^r B i( T )dT (4) 

0 0 

where Trf t P^(t) is the probability density function that the initial spin-down electron has relaxed to spin-up before 
tunneling out (recall we consider here the case that Prft = 0) = 0). 

The spin-down fidelity for dots with j > 1 are additionally affected by relaxation during the previous read-out 
stages including the emptying of all dots < j , and the time it takes to shuttle from their position in the array to 
the read-out position, which adds a shuttling time. The total waiting time before the j th dot is read-out is denoted 
by T° wait . This relaxation occurs with probability rf = J^ wait rf-e~ t ^ T ^dt. Therefore the spin down fidelity for the 
quantum dots with j > 1 is 

Fj — (S 7 + A J ')(1 — rf) + rf cr 7 (5) 

The last term accounts for the case that the spin down electron of the j th dot relaxes to spin up during T 3 wait but in 
the end is still detected as spin down. For the experiment the following values were used: T ^ = Trf= 130 /us, Tjj=300 
/us, T'.f valt = 222.15 /us and T^ ait = 514.3 /us. These values were chosen to optimize the sum of the three spin-down 
fidelities, whilst keeping the spin-up fidelity high. 

Note on the choice of T[: as we are transferring some of the spin states throughout the array, they will also 
experience different Pi’s depending on their location. In the current scheme each spin is waiting most of the time in 
its starting dot, so for the calculation of rf we take T(. During the read-out itself, so for the calculation of S - 7 and A - 7 
we always use T( = Trf Fig. ^l] shows the three measured T[. 


2.2 Determination of the tunnel rates 

For calculating the spin-down read-out fidelities, various tunneling rates need to be known. We can obtain T'rf t and 
r rf from the measurements of the histograms of the time it takes a spin-down electron to tunnel out (Tj ) and of the 
time it takes a spin-up electron to tunnel back into the quantum dot (T e 7 ), see Fig. ^2j We estimate T from the 
measured value of a- 7 , which characterizes the spin up fidelities and can be expressed as 


Ti- 0.7 fis 


rtd e - T lit 

1 out e 


dt 


r 


B- 7 (r) dr + i 


( 6 ) 
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Waiting time in (1,1,1) (ms) 


Fig. S 1: Single dot relaxation as a function of waiting time in (1,1,1)- Each datapoint is an average of 2000 
measurements. Grey lines are a fit to p 5 • e _t//T i + a?. 




Dot 3 


T d (|is) 


T d (|is) 


8000 


co 6000 I, r ’ = 75.8 kHz 
I' out 
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T e ^ 


T e ^ 
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Fig. S 2: Top row: histograms showing the distribution of the time Td it takes a spin-down electron to tunnel out at 

the read-out position of each respective dot for the measurement shown in Fig. The line is an exponential fit from 

which we determine the decay rate given by + -E. Bottom row: histograms showing the distribution of the time 

^1 

T e it takes a spin-up electron to tunnel back into the empty dot. The grey line is an exponential fit from which we 
can extract the decay rate given by F]^ J . 


Here fJ is the signal processing error which is found to < 0.1 % (to be discussed in the section on the characterization 
of the measurement bandwidth). 
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2.3 Characterization of the measurement bandwidth 


The signature of a spin-down electron tunneling out of the dot followed by a spin-up electron tunneling back into the 
dot is a pulse in the SD signal. Due to the finite measurement bandwidth, we can only detect pulses of sufficient 
duration and we therefore miss a fraction of the actual tunnel events. We characterize the measurement bandwidth 
by simulating such electron tunneling events by applying a rectangular pulse to P 3 with duration r and amplitude A. 
The triple dot is tuned deep into Coulomb blockade such that the applied pulse cannot result in a charge transition 
of any of the quantum dots. Due to the capacitive coupling of P 3 to the SD, a rectangular pulse will be created in 
the SD signal. The amplitude A is chosen such that the height of the pulse in the SD signal equals the signal from 
electron tunneling events for each respective read-out stage. The pulse duration is varied from 0.0 to 1.0 ps. Fig. £[3] 
shows the probability of detecting an event of duration r, P-'(r), as a function of the pulse duration for each of the 
three read-out channels. Using the formulas derived in Analytic expression for the fidelity we calculate the fidelities 
numerically by using the measured probabilities P - 3 (r); the result is shown in Table S |T] The fraction of the detected 
events at r = 0 s gives the signal processing error, e °, which corresponds to the fraction of traces in which an event is 
detected in the SD signal due to noise in the SD-signal during the read-out stage. 



0.3 0.4 0.5 0.6 0.7 

Pulse duration x (ps) 


Fig. S 3: The probability of detecting a pulse with duration r for each of the three read-out channels (P 3 (r)). Every 
read-out stage is monitored by a separate input channel of the FPGA. The output of the demodulation box is low-pass 
filtered at 1 MHz. The sensitivity of the SD to the three different read-out stages is slightly different as they occur at 
different detunings of the plunger gates. Each datapoint is an average of 1000 measurements. 


dot nr. 

Xi (ms) (worst,best) 

Spin-down fidelity (%) (worst,best) 

Spin-up fidelity (%) (worst,best) 

Orbital splitting (meV) 

1 

12.2 (11.4, 13.1) 

94.1 ( 93 . 8 , 94 . 5 ) 

97.4 (96.7, 98.0) 

2.0 

2 

11.5 (10.9, 12.2) 

95.8 (95.6, 96.1) 

99.3 (98.8, 99.9) 

1.8 

3 

8.5 (8.0, 9.0) 

97.4 (97.3, 97.5) 

98.0 (97.6, 98.4) 

1.0 


Table S 1: Read-out fidelities per dot. Values in brackets show the error margins based on the fits. For completeness 
the orbital splitting in each of the three dots was also measured using pulsed spectroscopy . 
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3 Detailed information of the applied pulse sequence 

In this section we give detailed information on the applied pulse sequences. Fig. ^4] shows the same charge stability 
diagram as Fig. lbc from the main text with additional labeling. Table S [2] gives an explanation including the details 
of the relevant stages. To correct for slow variations in the dot levels as a function of time (hours timescale), we always 
calibrate the three read-out stages before each longer measurement such as a complete T\ decay (~ 20.000 datapoints 
taken after one calibration run, which takes about half an hour). 



70 


50 


dV R F,SD | dV RFj SD 

dV R dV L 


(a.u.) 




(0,0,1) 

15 25 


Fig. S 4: Charge stability diagram for a M = -42 and b M = -56 from Fig. 1 from the main text with added details 
for the pulse sequence. 
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Stage(s) 

Purpose and details 

A 

Emptying stage (lasting 1 ms) during which all dots are emptied. 

B -4 C 

Loading an electron from the right reservoir into the right dot. We ramp from B to C in 600 /ns 
to load a '('-electron, or pulse from B to C to load an electron with random spin. 

C -4 D 

Shuttle electron from the right to the middle dot. Waiting time at D is 10 /ns. 

D -4 E 

Preparing to shuttle to the next dot. Waiting time at E is 1 /ns. 

E —>• F 

Shuttle electron from the middle to the left dot. Waiting time at F is 10 /ns. 

F -4 G 

Preparing to pulse into the other M-plane. Waiting time at G is 1 /ns. 

G -4 H 

Pulse into the other M- plane. Waiting time at H is 10 /ns. 

H -4 I 

Preparing to load a new electron into the right dot. Waiting time at I is 1 /ns. 

I -4 J 

Loading an electron from the right reservoir into the right dot. We ramp from I to J in 400 /ns 
to load a ('-electron, or pulse from I to J to load an electron with random spin. 

J -4 K 

Shuttle electron from the right to the middle dot. Waiting time at K is 10 /ns. 

K -4 L 

Preparing to load a new electron into the right dot. Waiting time at L is 1 /ns. 

L -4 M 

Loading an electron from the right reservoir into the right dot. We ramp from L to M in 400 /ns 
to load a ('-electron, or pulse from L to M to load an electron with random spin. 

M -4 N 

Optional stage to perform EDSR deeper into (1,1,1) to prevent photon-assisted tunneling with the 
reservoirs during the applied FM-burst. 

Read-out stage 3 

Gate voltage pulses larger than ~2 mV produce a spike in the RF-read-out signal. To prevent false 
events, we therefore first pulse to a position close to the read-out stage (2 mV more positive in R) 
and wait for 2 /ns. Only then, the RF-signal is unblanked and next we pulse into the read-out 
configuration for 130 /ns. 

L 

Empty the right dot for 70 /ns. 

L -4 K -4 J 

Shuttle the center electron to the right dot. Waiting time at J is 10 /ns. 

Read-out stage 2 

Similar as Read-out stage 3. 

I 

Empty the right dot for 100 /ns. 

I -4 C 

Path of I to C through all intermediate points with similar times as for the loading sequence. 

Read-out stage 1 

Similar as Read-out stage 3, except we now stay in the read-out configuration for 300 /ns. 

Compensation stage 

At the end of the pulse sequence we add a compensation stage of 10-45 ms that ensures that the 
total DC-component of the pulse is zero. This prevents unwanted offsets of the dot levels due to the 
bias tees. 


Table S 2: Detailed explanation of the applied pulse sequence as described by Fig. ^4] Unless noted otherwise, we 
always apply pulses from one point to the other. The total duration of this sequence varies between 3.3 and 54 ms, 
not including the compensation stage at the end. 
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4 Estimation of the error rate during shuttling 

In the main text we give estimates for the error rate during shuttling caused by (1) charge exchange with the reservoir 
and (2) random spin flips caused by the hyperfine interaction. In this subsection we will derive the models used to 
quantify these errors. 

4.1 (1) Charge exchange with the reservoir 

Fig. ^5] shows the pulse scheme of the multiple shuttling experiment in more detail. The probability that pulsing from 
A to B will be successful (i.e. the electron is fully transferred to B) can be calculated as follows: 

pt shuttle 

PA^B{t shuttle) = / T A^, B e~ TA ^ Bt dt (7) 

Jo 

where tshuttle is the time the electron waits in B before going to C, and is the interdot tunnel rate between 

the middle and right dot at that specific detuning. In this experiment the waiting time at each point A, B, C, D is 

equal tO tshuttle — 1 5-nhops * 



Fig. S 5: Pulse scheme for the multiple shuttling experiment. We first load a random spin at position A from the 
right reservoir. Next we go back and forth to A by applying the following pulse sequence: A—»B—»C—»D—)-C—»B—>A. 
Each run of going back and forth corresponds to four interdot tunnel events. The black dashed lines are extensions of 
the dot-reservoir charge transition lines. 

In the case that the electron was not successfully transferred to (0,1,0) at B, an error might occur when we continue 
pulsing to C. Pulsing from B to C crosses the extension of the right dot-reservoir charging line (dashed vertical line 
in Fig. # This implies that the following two processes can occur to the electron that is still in the right dot whilst 
the detuning has already been pulsed to C: (1) the electron will shuttle from the right to the middle dot with rate 
(2) the electron will tunnel out to the right reservoir with rate T reser voir, r, and a new electron is eventually 
loaded, but its spin will not be correlated with that of the initial electron, introducing an error. The probability that 
an error occurs going from B to C can therefore be expressed as: 

p _ f ^ reservoir,R 

* error,B—>C p -p 

J- A—>C i J- reservoir,R 

Leading to: 

Perror,A—>C — (1 — Pa—>b{ 1 shuttle)) ' Perror,B^>C 

In the following, we will neglect the possibility that none of the above two processes occurred at C. We will also 
neglect the possibility of an unsuccessful tunnel event from C -A D. This would lead to a second-order correction to 
the error in the reversed pathway D —>• B. The reversed pathway can then be expressed in a similar way as eq. ([9]). 
We measured that T reservo ir,R ~ 20 kHz, T reser voir,L ~ 10 kHz and Pd-^b > r> 1 MHz (lower bound, limited 
by the measurement bandwidth). The error should therefore be dominated by the path going from A to C. To get an 
upper bound for the error, we will assume from now on that P err or = P e rror,A^c = Perror,D^B- 

To calculate the measurement outcome at the end of the shuttling sequence, we only have to keep track of the last 
error in each measurement run. If an electron is exchanged for example twice during a run, we will eventually only 
measure the spin state of the last electron that entered. 


( 8 ) 

(9) 
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Depending on when the last error occurred, we load a spin-down electron with probability Pi oa d, 4 .- Its spin state will 
decay with Tf = Ti tWe i g hted- Please recall that P err or is equal to the error per two hops as we neglect errors in the 
processes C —t D and B —> A. 


Pl(n hops ,t C CD) = Pioad ,1 ((1 - P error ) nk °* s / 2 e- t ° CI >/ T ^ h ^') + 


k — Pt'hops /2 

Pioad,l- | ^ ^ Perror {1 Perror ) 

fc=l 


k-l e -(tcCD- nhov j 2 




( 10 ) 


The first term describes the process of no error occurring. The second term adds up all the possible errors starting 
from the last error happening in the final stage, k = 1, till the first stage, k = nhops/2- 


It was only possible to measure a lower bound for T^-^c > 1 MHz. Fig. Sjbjplots the outcome of eq. (10) for the 
same values of tccD as in the main text, Fig. 3a, for two different values of = F d^b'- 1 MHz and 2 MHz. Each 

curve is fitted linearly to extract a change in the spin-down fraction per hop. This linear approximation is only valid 
for our current regime of small error rates. For large error rates the spin-down fraction will increase non-linearly and 
eventually form a plateau at the value of Pi oa d,i- 

The fitted change in the spin-down fraction per hop is < 10 -8 for tccD > 7 ms, so we only quote the values for 
tccD = 2 ms: 18 • 10 -6 (T^^c = 1 MHz) and 0.5 • 10 -6 = 2 MHz). The experimentally measured values are 

given in Table S [3] 


tccD (ms) 

change in spin-down fraction per hop 

2 

-1.7 (-25,21) -10~ b 

7 

4.6 (-15,24) -10~ b 

35 

8.3 (-0.97,18) T0 _b 


Table S 3: The measured change in spin-down fraction per hop measured from the data in Fig. 3a) of the main text. 
Values in between brackets indicate the 95% confidence interval of the linear fit. 


Based on this simulation we can conclude that charge exchange with the reservoir is a negligible effect for tccD = 
7 and 35 ms and the probability to successfully shuttle from A to B or D to C will approach 1. For the case of tccD = 
2 ms, the slope extracted from the simulation is within the error bars of the experimentally measured value. It is not 
clear whether charge exchange with the reservoir is already the limiting mechanism here. Finally, we note that when 
shuttling electrons along longer arrays, charge exchange with the reservoirs will be negligible, if most dots in the array 
are in fact not coupled to reservoirs. 


a 


b 




nr. of hops 


nr. of hops 


Fig. S 6: Outcome of eq. 

a 1 MHz and b 2 MHz. In both cases r reservoir r = 20 kHz. 


(10) for the same values of tccD as the main text Fig. 3a for two different values of 

Each curve is linearly fitted (dashed 


lines) to extract a change in the spin-down fraction per hop. 
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4.2 (2) Random spin flips caused by the hyperfine interaction 


To estimate the error rate caused by the hyperfine interaction we describe what happens to a spin that starts in the 
spin-up state in the first dot. A similar reasoning applies to a spin that starts as spin-down. 

The spin-up state in the first dot can be described as a Bloch vector with length 1 and an angle 9 1 = arctan(Bjvi.j_/^ext) 


w.r.t to B ext where Sjvi,j_ = 



ss 5 mT is the magnitude of the local nuclear field orthogonal to B ext inside 


dot 1. In the next dot, the spin will experience a different nuclear field: B^ 2 ,±- The Bloch vector will now start 
rotating along a new axis given by 62 = arctan(B/v 2 ,j_/-®ea;t) w.r.t to B ext . In the current experiment, every step from 
one dot to the next takes longer than T|. Therefore the Bloch vector will soon dephase into a vector with length 
cos ($2 — 9i) and point at an angle of 9 2 w.r.t B ext . For small angles A 9 = 9 2 — 9\ ss arctan((i?jV2,_L — Bmi,a_)/ B ext ) = 
arctan((v^2-Bjvi,-L.)/-® ea;t ). Repeating this rihops times gives a final Bloch vector with length (cos(A 9)) nhopa at an angle 
^finaidot- The probability to be in the ground state of the last dot, and thus to be measured as spin-up, is then given 
by: 


P t = 0.5 + 0.5 • (cos(A 9)) nho *° 


1 Q _ (1 - cos(Afl)) • n hops 
2 


( 11 ) 


In the last step of the equation we apply the approximation that cos(Ad) r 3 1. The error per hop is then estimated 
by (1 - co ^ (Ae)) = 1 • IQ" 6 . 


5 Improvements to increase the fidelities of the current GaAs experi¬ 
ment 

This section covers the details of the described improvements for the current GaAs experiment as mentioned in the 
Discussion of the main text. These improvements combined give the red curve of Fig. 3c of the main text. 

To decrease the read-out time we can add a pulsed line to the SD. This allows us to always set the SD to the 
most sensitive point for each separate read-out position, which should easily improve the measurement bandwidth by 
a factor of two and thereby reduce the required read-out time by a factor of two. To decrease the emptying-tinre we 
propose to add a pulsed line to RS , so the tunnel rate between right dot and reservoir can be switched between fast 
(for emptying) and moderate (for read-out). Combined with larger interdot tunnel couplings both the emptying and 
the shuttling can take place on the ns timescale. Much faster shuttling would violate the adiabaticity condition, i.e. 
excitations to higher orbital states would occur [5J. The final improvement involves lowering the magnetic field to 
3.0 T (which would require to slightly reduce the electron temperature in order to maintain high read-out fidelities), 
which already doubles the T± time [5]. 
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6 Virtual gates L, M and R and the usage of ‘fast honeycombs’ 


‘Virtual gates’ 

To control the electron number inside each dot we change the voltages on gates Pi, P 2 and P 3 . In practice each gate 
couples capacitively to all three dots. Changing for example Pi which couples mostly to dot 1, will therefore also 
influence dot 2 and 3. To make selective control of each dot easier it is convenient to correct for this cross-capacitive 
coupling. This can be done by measuring the cross-capacitance matrix for the three gates recording the influence of 
each gate on every dot. Inverting this matrix allows you to create honeycomb diagrams with vertical and horizontal 
charge transitions in a so-called ‘virtual gate’-space. In such a ‘virtual gate’-space, the real gates Pi, P 2 and P 3 have 
been replaced by linear combinations of the three gates which allow the user to change the chemical potential in one 
dot, without changing it in the neighboring dot. 

For this experiment the exact relation between Pi, P2 and P3 and the ‘virtual gates’ L , M and R used is given by: 

/ L \ ( 1.56 0.0 0.40 \ ( Pi \ 

M = 0.0 1.56 0.0 P 2 

\ R J \ 0.19 0.51 1.05 J \ P 3 J 

This set of virtual gates did not perfectly correct for cross-capacitances, as the capacitances turn out to slightly 
vary themselves with the gate voltages. Working with virtual gates does however give one large freedom to take slices 
at any angle through the 3D honeycomb diagram that satisfy the requirements for the experiment. 


‘Fast honeycombs’ 

All charge stability diagrams shown in this work have been taken in so-called ‘fast honeycomb’ mode. Using the 
bias-tees connected to Pi, P 2 and P 3 it is possible to step one of them ‘slowly’ using a DAC and apply a triangular 
ramp on the other using the AWG. This significantly speeds up the measurements compared to stepping both gates 
using DAC’s. We have always plotted the reverse sweep, i.e. the voltage on horizontal axis is swept from positive to 
negative (with a rate of 100 mV/44 ms). The fading of the middle dot charging line can be explained through the 
small interdot tunnel couplings. Fig. .^JT] depicts the same stability diagram as in Fig. 1 of the main text over a larger 
scan range. 
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Fig. S 7: Charge stability diagram for a M = -42 and b M = -56 shown over a larger range than in Fig. 1. 
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7 EDSR spectra of each dot and spin-down initialization efficiencies 

Fig. ^8]shows EDSR spectra for each dot measured in a similar way as in [7], Despite similar (^-factors, time-variations 
of the local nuclear field still give rise to stable configurations where we can selectively address dot 2 and 3 in the 
(1,1,1 )- r egime. To perform such a measurement, the read-out sequence is started and then we manually tune the 
EDSR-frequency until we find a value that only addresses dot 3. We managed to find configurations stable for more 
than 30 minutes. We suspect there is some interplay between the electron and nuclear spins that leads to these stable 
points jB]. This requires further study that is outside the scope of this paper. For completeness Table S [4] gives an 
overview of the measured g -factors in each dot. 

Although we use adiabatic passage to create the spin-down states, the initialization efficiency is still limited on 
average to 76%. We attribute this non-perfect efficiency to two causes: (1) to achieve selective addressing of dot 2 and 
3 (which are closest in ^-factor) we on purpose slightly detune from the ideal spin-orbit mediated EDSR (SO-EDSR) 
resonance condition to prevent cross-talk. (2) At 3.51 T, the three hyperfine-mediated-EDSR resonance conditions are 
detuned by only 26 MHz to 46 MHz from the SO-EDSR condition [2]. Each resonance frequency by itself is capable 
of (partially) inverting the spin during a frequency chirp. The used FM-depth of 60 MHz could thus sometimes span 
several resonances that partially cancel each other. 



Fig. S 8: EDSR spectra for each dot. We load Ttt inside (1,1,1) and apply a frequency-modulated(FM)-burst 
to adiabatically invert the spin(s). Afterwards we read out the spin state. Each datapoint is an average of 300 
measurements. FM-depth = 60 MHz, burst time = 500 gs. Power applied to the gate is -29 dBm (excluding 
corrections for coax attenuation). The frequency is stepped from 21.35 GHz to 20.95 GHz 


dot nr. 

9i 

93 

1 

-0.4304 (0.4292, 0.4316) 

-0.0001853 (-0.0002596, -0.0001109) 

2 

-0.4337 (0.4307, 0.4368) 

-0.0001695 (-0.0003595, 0.0000205) 

3 

-0.4344 (0.4328, 0.4361) 

-0.0001762 (-0.0002814, -0.0000711) 


Table S 4: Measured g -factor in each dot based on fits of the form f res = ^hbB 0 f p^g 2i) of 

the main text. Values in brackets indicate the 95% confidence interval. 
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8 T\ measurements performed for different charge states in the honey¬ 
combs 

To further verify the proper operation of the single spin CCD we have also performed Tf measurements where we 
do not vary the waiting time in (1,1,1) as shown in Fig. 2 of the main text, but instead change the waiting time in 
different charge states. In Fig. ^9]we measure the T) decay of the spins by first loading three random spins in (1,1,1), 
and next vary the waiting time in (i) (1,0,1) or (ii) (1,0,0). In case (i) the right spin has already been read out before 
the waiting time, so its spin state will not change as a function of waiting time. The spin that was originally in the 
center dot is moved to the right dot before the waiting time, so it will now experience a Tf time that is more similar 
to T-l. Case (ii) shows no time-dependence for the center and right spin, and does show the Ti decay as expected for 
the left dot. 



Waiting time in (1,1,0) (ms) Waiting time in (1,0,0) (ms) 


Fig. S 9: a Ti decay in (1,0,1). Fits to the data give T-f = 14.3 ms, T 2 = 8.2 ms. The center spin now spends most 
of its time in dot 3, and therefore also relaxes with the same rate as T) 1 as measured in (1,1,1). b T) decay in (1,0,0). 
Fit to the data gives T? = 13.9 ms. 
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